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We consider various random models (directed polymer, ferromagnetic random Potts model, Ising 
- - - spin-glasses) in their disorder-dominated phases, where the free-energy cost F{L) of an excitation 

of length L present fluctuations that grow as a power-law AF(L) ^ L" with the so-called droplet 
, exponent uj > Q. We study the tails of the probability distribution n(a;) of the rescaled free-energy 

' cost X — ^^j^J'^ , which are governed by two exponents {rj-,r]+) defined by In n(a; ±00) ~ — |a;|''± . 

The aim of this paper is to establish simple relations between these tail exponents {rj-,rj^) and the 
■ droplet exponent uj. We first prove these relations for disordered models on diamond hierarchical 

lattices where exact renormalizations exist for the probability distribution Il{x). We then interpret 
these relations via an analysis of the measure of the rare disorder configurations governing the tails. 
Our conclusion is that these relations, when expressed in terms of the dimensions of the bulk and 
of the excitation surface are actually valid for general lattices. 
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INTRODUCTION 



To understand the low-temperature phase of disordered systems, it is important to characterize both the statistics 
of the ground state energy over the disordered samples, and the statistics of excitations above the ground state within 
one sample. 

The ground-state energy Eq of a disordered sample is the minimal energy among the energies of all possible 
configurations. The study of its distribution thus belongs to the field of extreme value statistics. Whereas the case of 
independent random variables is well classified in three universality classes [l| , the problem for the correlated energies 
(~| within a disordered sample remains open and has been the subject of many recent studies. The interest lies both 
Q [ (i) in the scaling behavior of the average E§'"{L) and the standard deviation AEo{L) with the linear size L 
CJ . (ii) in the asymptotic distribution P{x) of the rescaled variable x — (i?o — Eq" [L)) / /S.Eq{L) in the limit L ^ 00 

00 : 

0^ ■ For spin-glasses in finite dimension d, a sample of linear size L contains N — L'^ spins. Following the definitions of 
' Ref. [3], the 'shift exponent' cug governs the correction to extensivity of the averaged value 

(N 

ES^iL)^L''eo+L^^ei + ... (2) 

f — ' Within the droplet theory [1,0], this shift exponent ujg coincides with the domain wall exponent ujdw and with the 
. droplet exponent lu of low energy excitations (see below) 

' OJs = i^DW = W (3) 

V : . 

■ The 'fiuctuation exponent' ujf governs the growth of the standard deviation 

l^Eo{L) - L'^f (4) 

In any finite dimension d, it has been proven that the fluctuation exponent \s ujf = d/2 [51. Accordingly, the rescaled 
distribution P{x) of Eq. ([T]) was numerically found to be Gaussian in d = 2 and d = 3 [4], suggesting some Central 
Limit theorem. In contrast with finite-dimensional spin-glasses where one needs to introduce two exponents and uo / , 
the directed polymer model @ is characterized by a single exponent uj that governs both the correction to extensivity 
of the average E^'"{L) and the width AE(){L) 

E^^iL) ^ Lea + L'^ei + ... 

AEn{L) -L^e2 + ... (5) 

This exponent also governs the statistics of low excitations within the droplet theory [3| , as confirmed numerically [3] ■ 
In dimension 1 -f 1, this exponent is exactly known to he Lj{d = 1) = 1/3 [1, H, [13, EI and the corresponding rescaled 
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distribution n(a;) is related to Tracy- Widom distributions of the largest eigenvalue of random matrices ensembles 

Among disordered systems, the directed polymer model thus presents the following distinctive feature : the statistics 
of the ground state energy over the samples is directly related to the statistics of excitations within a given sample 
In finite-dimensional spin systems however, the statistics of the ground-state energy over the samples is not very 
interesting : the fluctuation exponent 6f = d/2 and the corresponding Gaussian distribution P{x) simply reflects the 
fluctuations of the L'^ random couplings defining the samples. The only information it contains on the statistics of 
excitations is the shift exponent lUs = uj governing the correction to extensivity of the averaged value. In this paper, 
we will be interested into the statistics of the energy El of excitations above the ground state in a given sample. Its 
fluctuations are governed by the droplet exponent u 

AEl - (6) 

and its distribution is expected to follow a scaling form as L — > oo 

Our aim is to show that the exponents governing the tails of the rescaled distributions n(a;) are simply related to 
the droplet exponent uj. Since the whole low-temperature phase (T < Tc) is described at large scale by the zero- 
temperature fixed point, the statistics of the free-energy Fl of excitations at T < Tc is the same as the statistics of 
the energy of excitations above the ground state, up to some rescaling with the correlation length ^(T). So the 
results concerning the tails of Tl{x) also concerns the tails of the free-energy of excitations at any temperature T < T^,. 

To establish the relations existing between the tails of the probability distribution 11 (x) and the droplet exponent 
u, we will first focus on the diamond hierarchical lattices where exact renormalizations exist as we now recall. Among 
real-space renormalization procedures , Migdal-Kadanoff block renormalizations Ti*] play a special role because 
they can be considered in two ways, either as approximate renormalization procedures on hypercubic lattices, or as 
exact renormalization procedures on certain hierarchical lattices [H, [13] • One of the most studied hierarchical lattice 
is the diamond lattice which is constructed recursively from a single link called here generation n = (see Figure [IJ: 
generation n = 1 consists of b branches, each branch containing 2 bonds in series ; generation n = 2 is obtained by 
applying the same transformation to each bond of the generation n — 1. At generation n, the length L„ between the 
two extreme sites A and -B is L„ = 2", and the total number of bonds is 

B„ = {2br - Lt-'-^'^ with 4//(6) = (8) 
where deff(b) represents some effective dimensionality. 




FIG. 1: Hierarchical construction of the diamond lattice of branching ratio b. 



On this diamond lattice, various disordered models have been studied, such as for instance the diluted Ising model 
17 1 jferrornagnetic random Potts model [3 [l^, [2^ , spin-glasses [2l|, [U, [2^ and the directed polymer model 

26 . 127I [28l. [29l. [30I [Sll. [3^ [ssl. [33| . In this article, we start from the exact renormalizations existing for these disordered 
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models on the diamond lattices to derive the relations existing between the tails exponents and the droplet exponent 

LO. 

The paper is organized as follows. The tails of the distribution I\-{x) of the ground state energy of the directed 
polymer are discussed in Section together with numerical results; the tails of the distribution n(a;) of excitations 
in the ferromagnetic random Potts model and in the Ising spin-glasses are studied in Section IIIII and in Section IIVI 
respectively. Finally in Section |Vl we generalize these results to other lattices via an analysis of the measure of the 
rare disorder configurations governing the tails. Our conclusions are summarized in Section IVTl Appendix VK\ contains 
more detailed calculations, and Appendix [B] contains a brief reminder of Zhang argument for the directed polymer on 
hypercubic lattices. 



II. DIRECTED POLYMER ON DIAMOND LATTICE 



In this Section, we study the tails of the rescaled probability distribution 11 (x) for the ground state energy of the 
directed polymer model. (Eq. [T]). 



A. Reminder on the exact renormalization 



1. Renormalization at finite temperature 

The model of a direct ed poly mer in a random medium [6| can be defined on diamond hierarchical lattice with h 
branches [H, HI, [IE HE IS [HI, [H, [H, Is^l ■ The partition function Z„ of generation n satisfies the exact recursion 

M 

Z„+i=^Z(2-i)Z^) (9) 

a=l 

where {Zn \ zl^^^) are (26) independent partition functions of generation n. At generation n = 0, the lattice 
reduces to a single bond with a random energy e drown from some distribution p{e) and thus the initial condition for 
the recursion of Eq. [5] is simply 

Z„=o = e-^' (10) 

In the low-temperature phase where the free-energy width AF(L) grows with the scale L, the recursion is dominated 
at large scale by the maximal term in Eq. [9] 

Z„+i ~ maxi<a<b (11) 

or equivalently in terms of free-energies Fn = —TlnZn 

F„+i ~ fmni<a<b (^^P""'^ + F^'"') (12) 

This effective low-temperature recursion coincides with the recursion of the energy Eq of the ground state studied in 
[26L l27l|. The whole low-temperature phase is thus described by the zero-temperature fixed-point. 



2. Renormalization at zero temperature 



We now focus on the statistics of the ground state energy of the directed polymer [2^ |28|, [291, \3^. At T = 0, the 
recursion for the ground state energy involves the following minimisation [26| 



En+i ~ niin 



This translates into the following recursion for the probability Pn(E) [20| : 



/•+00 r I' + ca 

/ dEPn+liE)^ / dSQniS) 
J z lJ z 



(13) 



(14) 
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where QniS) is the distribution of the sum 

/ + 00 /• + 00 

dEi / dE2Pn{El)Pn{E2)5{S - {E^ + E2)) 
-00 J —00 



(15) 



3. Renormalization in the scaling regime 
For large n, one expects the scahng [2^ 

n^oo dn V On 

where the term 7„ is extensive in the length L„ 

7„ ~ eoLn 

71 — ^OO 



i.e. 7n+i/7n 2. The width 5„ scales with the length L„ with some a priori unknown exponent w [2£ 

(5„ ~ A" ~ with A = 2" 
Replacing the scaling form of Eq. \W\ in the recursion of Eqs [T4j [15] yields 



+00 



-Tti + 1 



/+00 r+oc 
dx / dxi / da:2nb(xi)nb(a;2)(^(a; - (xi + 0:2)) 
J —00 J —00 



Using 7„+i/7„ = 2 and (5„_|_i/(5„ ^ A = 2", one obtains 

f +00 

c?xnh(a;) = 



dxGb{x) 



where 



Gbix) 



+00 



dxi 



dx2llb{xi)Ilb{x2)6{x - {xi + X2)) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



The recursion simplifies in the limit 6=1, where it reduces to the Central Limit theorem for the sum of random 
variables with [291 



w(6=l) = 

n,,=i(x) = 



1 



(22) 



We refer to [26| where an expansion in = 1 + e has been developed. 

Another limit where the recursion simplifies is 5 — *■ 00. In the limit, a single iteration consists in taking the minimum 
of a large number 6 — > c» of random variables. The rescaled distribution is then the Gumbel distribution [IJ 



uj^b — + 00) = 



(23) 



For 1 < 6 < +00, the exponent uj{b) is expected to decay from uj{b = 1) = 1/2 to uj{b = 00) = 0, Accordingly, the 
rescaled distribution Ilb{x) is expected to interpolate between the Gaussian and the Gumbel distribution. 
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B. Relations between tails exponents and the droplet exponent 

Let us now focus on the tail exponents 77(6) and r]'{b) of the probability distribution 11^ 

lnnb(a;) ~ -c(-x)''('') + ... 

x—^ — oo 

\nUbix) ~ -dx'''^'') +... 



(24) 



From the two extreme cases of Eqs [52] and [531 one expects that the left exponent 77(6) varies between 

,7(6=1) =2 

7/(6 00) =1 (25) 

whereas the right exponent rj'ib) varies between 

ry'(6=l) =2 

r]'{b = 00) =00 (26) 
The aim of this section is to show that these tail exponents are simply related to the droplet exponent uj{b) via 

1 



ri'ih) 



1 - Lu{b) 

ln(26) 



In 2 



,y(6) 



(27) 



To make things clearer, we have chosen to present here in the text only simple saddle-point arguments at leading 
order. We refer to Appendix [X] (see Eqs IA14 1 and lABlj) for a much more detailed proof with subleading corrections. 

1. Left-tail exponent r] 

Assume that the left-tail decay of the probability distribution I\-b{x) is given at leading order by Eq. [24l A 
saddle-point analysis shows that the probability distribution fEg I2ip of the sum x — xi -\- X2 presents the tail 



l^Gbix) 

The recursion of Eq. [501 yields by differentiation that the tail of lib is related to the tail of Gf, via 



1 „ fu\ 
A^^ (a) u 



bGb(u) 



+ 00 



dyGbiv) 



6-1 



bGb{u) cx be-'^<-'iy+- 



Using A = 2*^, this yields in terms of the variable x = u/X 



\nUb{x) ~ _2cf-4^V + ... = -2i+"("-i)c(-x)'' + 



(28) 



(29) 



(30) 



The consistency with the scaling form of the tail of Eq[24| yields ry = 1/(1 — w) as stated in Eq [27l 



2. Right-tail exponent r\ 



Assume that the right-tail decay of the probability distribution n(,(a;) is given at leading order by (Eq. [24]). A 
saddle-point analysis shows that the probability distribution Gh (Eq ISTll of the sum x = xi + X2 presents the tail 



lnGfc(a;) 



- -2c f 

X — > + oo V 2 / 



The recursion of Eq. [20] yields by differentiation that the tail of Ilh is related to the tail of Gb via 



a"^ (a) 



bGbiu) 



+ 00 



dyGbiy) 



b-l 



(X [Gbiu)Y 



oc e 

u — ' — 00 



-2bc{-^y+... 



(31) 



(32) 
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Using A = 2", this yields in terms of the variable x = u/X 



(33) 



The consistency with the scaling form of the tail of Eq[M]yields 7y'(l ~ oj) = ln(26)/ ln2 as stated in Eg [771 



Discussion 



The relations of Eq. [^T] have already been found for the special case 6 = 2 in [2^ . However in [23| , a third relation 
between the three exponents was also written, leading to the simple values 77(6 = 2) — 3/2 and ri'{b = 2) = 3 that 
seem now excluded numerically (see below). 



C. Numerical results of the ground state energy statistics 

1. Method : numerical recursion of the probability distribution 



In AE„ 



2 




0.5 



Q)(b) 



0.4 



0.3 - 



0.2 




0.1 - (b) 



10 



InL 



20 30 
b 



FIG. 2: (Color online) Width AEo{L) of the ground state energy distribution for a polymer of length L (a) the log- log plot of 
width A^o ~ L'^^''^ for b = 2,3, 4, 5, 6, 7, 8, 12, 16, 24. (b) Exponent Lu(b) as a function of the branching ratio b (Precise numbers 
are given in Table |T| 

As explained in [261 , it is more convenient numerically to consider the iteration of discrete probability distributions 
of the form 



+00 



Pn{e) = ^ a„(m)5e 



(34) 



m— — 00 



where the energy e can take only integer values. This form is conserved via the recursion of Eq[T3]that only involves 
summation of energies and choice of minimal value. The convolution step of Eq. [15] can be written as 



Q„(e) = (P„ * P„)(e) = J2 bnim)Se 



(35) 



with the following rules 



-t-CXD 



bn{rn) ^ ^ an{m')an{m - m') 



(36) 



m — — 00 
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Since the function 



/• + 00 +°P 

/ deQn{e) = bn{m) = Cn{m)6{m — 1 < x < m) 



m— — oo 



is constant on intervals with values 



it is easy to raise it to power b 



+00 



+00 



Cn(m) ^ ^ bn{m) 



deQn{e) 



+00 



= [cn(m)]^ 0{Tn — 1 < X < m) 



Since we have 



m— — 00 



+00 



+00 



/ deP„+i(e) = ^ a„+i(m) = ^ ^ a„+i(m') 

^ m>2; m' = — 00 Lm'=m 

the recursion of Eq. [14] yields 

+00 

a„+i(m') = [c„(to)]'' 



0(m — 1 < a; < m) 



and thus the coefficients at generation [n + 1) can be obtained via 

an+i{m) = [c„(to)]'' - [c„(m + 1)]^ 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 



This method allows to obtain very accurate results for the fluctuation exponent Lu{b) and for the tail exponents 
r]{b) and r/(&) because the probability distribution can be evaluated very far in the tails. The results presented below 
have been obtained by the iteration up to Ug ~ 30 generations of the discrete distribution of Eq. [53] with numerical 
tail cut-offs of order —10000 < rn — mcenterin-) < 1000 around the moving center mcenter{n) for each generation n. 



2. Numerical results for the fluctuation exponent uj{h) 



b 


uoib) 


riib) 


2 


0.299 


1.43 


3 


0.236 


1.31 


4 


0.205 


1.26 


5 


0.186 


1.23 


6 


0.173 


1.21 


7 


0.163 


1.20 


8 


0.156 


1.18 


12 


0.141 


1.16 


16 


0.131 


1.15 


24 


0.123 


1.14 



TABLE I: Results for the exponents Lo{h) and rjib) as the number b of branches varies 

We first show on Fig. [5] the log-log plot of width Ai?o ^ L'^^'') of the ground state energy probability distribution. 
The measures of the slopes yield the values given in Table |TJ The values of w(6) are in agreement with the existing 
previous numerical measures [H, [H, |3l|- The curve w(&) shown on Fig. seems to suggest that the droplet 
exponent u) remains positive as long as the effective dimension deff{b) — ln(2&)/ln2 (Eq[8]) remains finite. Note 
that for hypercubic lattices, the existence of a finite upper critical dimension dc above which the droplet exponent 
vanishes has remained a very controversial issue between the numerical studies [35l |3^ [37l . [38j and various theoretical 
approaches [H, SO, |4l| . 
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FIG. 4: (Color online) Log-log plot of the tails of the rescaled probability distribution of the ground state energy for branching 
ratio b = 2,3,4,5,6,7,8,12,16,24. (a) Measure of the exponent 77(6) governing the left tail lnlli,(a; — > —00) ~ —{—x)^^''^ (b) 
Measure of the exponent ri'{b) governing the right tail In Ilb{x — + +00) ^ '''' The curvatures of ln(— In lib (x)) as a function 
of ln(x) show that the corrections to the leading behavior are stronger than for the left tail shown in (a). 



3. Numerical results for the rescaled probability distribution fit (a;) 



We show on Fig. [3]the asymptotic rescaled probability distributions Ilb{x) defined in Eq. [T] On Fig. [3] (a) we show 
Ilb{x) in the bulk for b — 2,A, 8. On Fig. [3] (b) we show lnnh(a:;) to see the behaviors far in the tails. 

From the Figures HI we have measured the tails exponents defined by Eq. [24l The left-tail exponent ri{b) can be 
measured as the slopes of the curves of Fig. |l](a). The results 77(6) are given in Table[Tl The relation 77(6) = 1/(1— uj{b)) 
(Eq. [27]) is well satisfied by our numerical results. 

The right-tail exponent 77' (6) turns out to be more difhcult to measure precisely, because the corrections to the 
leading behavior are much stronger, as shown by the curvatures of our results of Fig. [4] (b). However, the values 
predicted by the relation of Eq. [57] 

V'{b) = 77(6) (43) 



are compatible with our data. 
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D. Extension to the critical point 



On the diamond hierarchical lattice, the free-energy fluctuations of the directed polymer either grow as for 
T < Tc, or decay as for T > or remain of order 0(1) exactly at T^. The study of the tails of the critical 

rescaled distribution [i^ is a special case of the above relations fEq [77]) with lUc — 

Vcib) = 1 

-'^ m 

We refer to Eal for more details. 



III. ISING AND POTTS RANDOM FERROMAGNETS ON DIAMOND LATTICE 



In this section, we study the tails of the rescaled distribution n(a;) for the excitations in random Ising and Potts 
ferromagnets (Eq. [7]). 



A. Reminder on the exact renormahzation for the effective coupling 

The Ising Hamiltonian reads 

H J sing — ^ ^ '^ij^i^j (45) 

<i,3> 

where the spins take the values Si = ±1 and where the couplings J^.j are positive random variables (see below section 
IIVI on spin-glass for couplings of random sign). The Potts Hamiltonian is a generalization where the variable ai can 
take q different values. 

Hpotts - - ^ 2J,J,^^„^ (46) 

<h3> 

(We choose (2 J) to recover Ising for q = 2) 



1. Renormalization at finite temperature 



The effective coupling between the two end-points A and B of the diamond lattice of Fig. [T] is defined by 



Z++ 



(47) 



where ^++ and Z_| are the partitions functions corresponding respectively to the same color at both ends or to two 

different colors at both ends. So the effective coupling represents the free-energy cost of creating an interface between 
the two ends at distance L. The renormalization equation reads in terms of the variable y = e^^'' jl8i, |19|, S3| 



Vn+l 



y^'^ + {q-'i), 



(48) 



In the low-temperature phase where the effective couplings J„ grow with the length L„ = 2", the contribution of each 
branch (i) is dominated by the minimal coupling between {Jn^\ Jn^^)- This leads to the effective zero-temperature 
recursion 



J n 



+ 1 



7(ii) j{^2) 



(49) 



The whole low-temperature phase is thus described by the zero-temperature fixed-point. 
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2. Renormalization at zero temperature 



We now focus on the zero-temperature renormalization 

b 

J„+i=^min [4'^\J(*^) 



i=l 



(50) 



Note that two operations 'sum' and 'min' occur in the opposite order with respect to the directed polymer case of 
Eq. [T31 Eq. [501 translates into the following recursion for the probability P„(J) : 

p + oc r+oc 

Pn+iiJ) = / dK,Pr''{K,)... / dK„Pr'^{Kb)5{J - {K, +K2 + ... + Kb)) (51) 
Jo Jo 

where P™™(_ftr) is the distribution of the minimum K = min{Jn\ j'^^) of two variables drawn with the distribution 

Pn{J) 

2 

(52) 



p-\-oo r-\-oo 

/ dKPr''iK)= / djp„(j) 

J z z 



3. Renormalization in the scaling regime 



For large n, one expects the scaling 



Replacing the scaling form of Eq. [53] in the recursion of Eqs [5T1 [52] yields 



dxinr"(a;i). 



dxbll^'''{xb)6{u ~ (xi + X2 



and 



da;n™"(a;) 



+ 00 



dxUbix) 



-Xb)) 



(53) 



(54) 



(55) 



and where the term 7„ should grow asymptotically as jn+i/ln ^ b. In terms of the effective dimension df,ff{h) of Eq 
[U this corresponds to 



d,ff(b)-i 



(56) 



as it should for an interface of dimension dg = deff{b) — 1. 

The recursion simplifies in the limit b —^ oo ^2'3,\ where a single iteration consists in summing a large number of 
random variables 

1 _x2 

nh_.oo(a;) = ^=e 2 (57) 



/27r 

(Note that here, the usual other simple limit 6 = 1 is not interesting since it does not correspond to a positive droplet 
exponent). 



B. Relation between the tail exponents and the droplet exponent 



Let us now focus on the tail exponents 77(6) and ri'{b) of the probability distribution Hb defined as in Eq[24] The 
aim of this section is to show that these tail exponents are simply related to the droplet exponent w(6) via 



Vib) 



V'{b) 



In 2 



ln(2b) 
In 2 



(58) 



As in the directed polymer case, these relations can be obtained via simple saddle-point arguments. 
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1. Left-tail exponent rj 

If the left tail of lib (a;) is described by Eq[24l the distribution n™«" presents the same decay by differentiation of 
Eq. [m 

n™"(a;) = 2nb(a;) / dyUb{y) ~ 2nb(x) ^ 2e-<--^'' (59) 



A saddle-point analysis shows that the convolution of b variables Ki then presents the following decay (Eq [ 

in.(^) . e-K-f)" (60) 
Using A = 2'^, this yields in terms of the variable x — u/X 

InUbix) ~ ^bc(-^y ~ -6i-"2""c(-a;)'' (61) 



a: — > — oo 



The consistency with the tail of Eq [24] yields the following constraint 2^'^ = 5'' ^ leading to the result for ri{b) given 
in EqESl 

2. Right-tail exponent rj 

If the right tail of 11^ [x) is described by Eq [24l the distribution 11™*" presents the following exponential decay by 
differentiation of Eq. 

nr"(a;) -2nb(a;) / d2/nb(y) cx H^Ox) cx e^^rf^" (62) 



A saddle-point analysis shows that the convolution of h variables Ki then presents the following decay (Eq [54 

1~ /^^\ -2^-)"' 



-nb(^^) c._e— "^l*^ (63) 



Using A = 2'^, this yields in terms of the variable x = u/\ 

/Ar\''' , , , 

lnnh(a;) ~ -2hd\ — \ ~ -51-" 21+" "da;" (64) 

The consistency with the tail of Eq[M] yields the following constraint 2^+'' ^ = W leading to the result for 77' (6) 
given in Eg 1581 

C. Extension to the critical point 

As in the directed polymer case (Eq [33]) , one may try to extend the results on the tail behaviors in the low- 
temperature phase where w > to the critical point where = 0. The right tail exponent becomes at criticality (Eq 
[51]) 

Vcib) = (65) 

For the left tail however, the problem is qualitatively different at criticality and will not be discussed here ( the critical 
invariant distribution Pc{J) does not extend to (—00) anymore, but reaches the natural boundary J = 0.) 

IV. SPIN-GLASSES ON DIAMOND LATTICE 

We now consider the case of Ising spin-glasses described by Eq. [JS] but now the random couplings J^j can be 
positive or negative and are distributed with a symmetric distribution Po{J) 

Po{J)^Poi~J) (66) 
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A. Reminder on the exact renormalization for the effective coupling 

1. Renormalization at zero temperature 

The renormalization at finite temperature is still described by Eq. [48] with q = 2. However, the fact that the 
couplings are of arbitrary sign yields that the zero-temperature renormalization now reads 



J„+i = ^ sign (j^) J('^)) min | J^) 



(67) 



instead of Eq [501 corresponding to the case of positive couplings only. The symmetry of the initial condition (Eq Wf} 
is conserved by the renormalization 



Pn{J)=Pn{-J) 



(68) 



To translate the renormalization of Eq. [57] into a recursion for the probability distribution Pn{J), it is convenient to 
introduce the probability distribution P"'"' of | J„ | 



PriJ >0)^ Pn{J) + P„(- J) = 2P„(J) 



and the probability distribution Q'^^{K) oi K = min 



> satisfying 



+ 00 



dKQT{K) 



+ 00 



djpriJ) 



(69) 



(70) 



(the value z — Q corresponds to the normalization condition of both distribution). The probability distribution Qn{J) 



oi K — sign (^Jn^Jn ^^ mil 



> is symmetric in K —K 



Qn{K ^ 0) = -gf ^ {\K\ 



The recursion of Eq. [S7| corresponds to the convolution 

/+CX; /•+00 
dKiQn{Ki)... / dKbQn{Kb)5{J - (Ki +K2 + ... + Kb)) 
-00 —00 



(71) 



(72) 



2. Renormalization in the scaling regime 



For large n, one expects the scaling 



p(j) ~ J-Hb ( — 



(73) 



where 7„ = in contrast with the ferromagnetic case of Eq[53l and where the scaling function 11^ is symmetric ( Eq. 



Ilbix) =Ilb{-x) 



(74) 



B. Relation between the tail exponent and the droplet exponent 

Note that the true rescaled distribution 11*'""'^ (x) as defined by Eq. [Tjwith a positive energy for excitations above 
the ground state actually corresponds to the distribution of the absolute value of the rescaled coupling lib (a;) defined 
in Eq[73| 



n'r'{x) = nb{x) + nbi-x) 



(75) 
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So the true distribution Il*J^'^{x) begins at x = with a non-zero value n*™'=(0) H and there exists a single tail 
as a: ^ +00. So in contrast with the previous cases of the directed polymer and of the ferromagnetic Potts model 
described by two tail exponents (Eq[24|), the distribution Ui,{x) of the rescaled coupHng for spin-glasses presents a 
single exponent rj'ib) as a consequence of the symmetry of Eq. [71] 

Inllbix) ~ -d\x\'^' (76) 

X — ^±00 

We will not repeat here the calculations that are similar to the case of the right tail of the Potts model and that lead 
to the same relation between exponents fEg [55)1 

ln(2b) 

-^w 



V. GENERALIZATION OF THE TAIL EXPONENTS RELATIONS TO OTHER LATTICES 

In the previous sections, we have derived the relations that exist between the tail exponents {ri,r]') and the scaling 
exponent u for various disordered models on the diamond hierarchical lattices from the exact renormalization recur- 
sions. The obtained relations are actually very simple in terms of the effective dimension deff(b) of these lattices (Eq. 

This suggest that these relations should have a simple interpretation. In this section, we explain the physical 
meaning of these relations and generalize them to other lattices. 



A. Tail exponents for the directed polymer 

1. Physical interpretation of the left tail 

The relation rj = 1/(1 — to) { Eq. [27]) derived previously from the exact renormalization on the diamond lattice 
can be interpreted as follows. The left tail of the ground state energy of the directed polymer corresponds to samples 
that leads to much lower energy than the average. Let us evaluate the probability to obtain a ground state energy 
E/Q — 7n — o,Ln extensively below the averaged value 7„ of the scaling function of Eq. [Tn]with tail behavior described 
by EqM 

Pn{Eo=ln-aL^)c^^nJ-a^) cc e-^-^"^"""' (78) 

\ KJ ^^^^ 

On the other hand, to obtain such a ground state energy Eq — 7„ — aL„ extensively below the averaged value 7„, it 
seems reasonable to ask that each bond of the ground state configuration of length L should have an energy ti lower 
than the average, which happens with the exponentially small probability 

L 

~ e"""*-^ (79) 

The identification of the length exponents in Eqs [78] and [79] yields 

r]{\-uj) = l (80) 

which corresponds to the relation found for the diamond lattice ( Eq. I27p . From this interpretation in terms of the 
rare disordered samples that govern the left tail of the ground state energy configuration, we expect that the relation 
of Eq[80]is actually valid on any lattice, and in particular for hypercubic lattices oi 1 + d dimensions. The relation 
is satisfied by the exact exponents in 1 + 1 dimensions with lo = 1/3 [1, [M El and 77 = 3/2 0, [H [II- The 
relation of Eq[5n]has been previously derived for hypercubic lattices via Zhang argument ^ (the argument is recalled 
in Appendix [B] for comparison) and has been checked numerically in [43(| for d = 2, 3. However Zhang argument only 
concerns the left tail because it is based on the existence of a Lyapunov exponent for positive moments of the partition 
function. It cannot be extended easily to the right tail that would be in correspondence with negative moments. This 
is in contrast with the rare events analysis presented here that can be extended to the right tail as we now explain. 




dep{e) 
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2. Physical interpretation of the right tail 



The relation 77' = def f{b) / {1 — lj) ( Eq. [27]) derived previously from the exact renormalization on the diamond 
lattice can be interpreted as follows. The right tail of the ground state energy of the directed polymer corresponds to 
samples that leads to much higher energy than the average. Let us evaluate the probability to obtain a ground state 
energy = + aLn extensively higher the averaged value 7„ of the scaling function of Eq. [16] with tail behavior 
described by EqIMI 



oc 



(81) 



On the other hand, to obtain such a ground state energy -Bo = 7n + o-Ln extensively higher the averaged value 7„, it 
seems reasonable to ask that all L'^''ff bonds of the lattice should have an energy higher than the average, which 
happens with the exponentially small probability 



dep{e) 



-cstL'^^ff 



The identification of the length exponents in Eqs [HU and [55] 

Tj'{l - Cj) = deff 



(82) 



(83) 



exactly corresponds to the relation found for the diamond lattice ( Eq. [77]) . From this interpretation in terms of the 
rare disordered samples, we expect that the relation of Eq [83] is actually valid on any lattice, and in particular for 
hypercubic lattices of dpjf = I + d dimensions. The relation is satisfied by the exact exponents in 1 + 1 dimensions 
with uj = 1/3 [1, i, E [ill and r/ = 3 [13, [Hill- 



B. Tail exponents for the ferromagnetic random Potts model 

1. Physical interpretation of the left tail 



The relation derived for the left tail of the Potts model on the diamond lattice fEa[55)) reads in terms of the effective 
dimension dej / (6) of Eq [5] 



Vib) 



deff{b)-l 

deff{b)-l-^{b) 



(84) 



We now propose the following physical interpretation. The left tail corresponds to samples that leads to much lower 
effective coupling than the average 7„ ~ L^^"^^''^' ^ (Eq. Let us evaluate the probability to obtain an effective 



coupling J — — aL 
behavior described by Eg [24] 



extensively below the averaged value 7„ of the scaling function of Eq. [53] with tail 



P„(J = 7n - aLt^"^'^-') ^ -i^n, -a "^" 



oc e 



(85) 



On the other hand, to obtain such a low effective coupling J = jn — aL!^^^^^^ ^ extensively below the averaged value 
7„, it seems reasonable to ask that each bond of the interface of dimension should have a coupling Ji lower 

than the average, which happens with the exponentially small probability 



dJPo{J) 



stL'^'lf 



(b)~i 



The identification of the length exponents in Eqs [55] and [55] 

V{deff{b) - 1 - w) = deff{b) - 1 
exactly corresponds to the relation found for the diamond lattice ( Eq 



(86) 



(87) 
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2. Physical interpretation of the right tail 

The relation derived for the left tail of the Potts model on the diamond lattice fEglSS]) reads in terms of the effective 
dimension deff{b) of Eq[5] 

deff[b) - l-uj[b) 

We now propose the following physical interpretation. The right tail corresponds to samples that leads to much higher 

effective coupling than the average 7„ ~ hft'^^^'^ ^ (Eq. [56]) Let us evaluate the probability to obtain an effective 

coupling J = 7n + aLn'"'"^'' ^ extensively above the averaged value 7„ of the scaling function of Eq. [53] with tail 
behavior described by Eq [24l 

PniJ = 7n + aLt^"^''^-')^i-Il,\a^^, I oc e-''''"^" (89) 




On the other hand, to obtain such a high effective coupling J = 7„ + aifi"^^^''^ ^ extensively above the averaged value 
7„, it seems reasonable to ask that all bonds of the sample of dimension L'^^/zC') should have a coupling Ji higher 
than the average, which happens with the exponentially small probability 



dJPoiJ) 

The identification of the length exponents in Eqs [55] and [TO] 



f]'{d,ff{b)~l-iu)^d,ff{b) (91) 

exactly corresponds to the relation found for the diamond lattice ( Eq. 



VI. CONCLUSION 

In this paper, we have studied the statistics of excitations of finite-dimensional random models (directed polymer, 
ferromagnetic random Potts model, Ising spin-glasses) in their low-temperature phase characterized by a positive 
droplet exponent w > 0. We have shown that the tails of the rescaled probability distribution are characterized by 
two tails exponents {r],T]') that are simply related to the droplet exponent uj. We have first proved these relations on 
the diamond hierarchical lattices where exact renormalizations exist for the rescaled probability distribution. We have 
then given the physical meaning of these relations in terms of the measure of the rare disorder configurations governing 
the tails. This interpretation allows to understand the asymmetry rj < r]' because a 'good' sample contributing to 
the left tail is a sample containing L"^' 'good' random variables for an interface of dimension dg < d, whereas a 'bad' 
sample contributing to the right tail is a sample containing L'^ 'bad' random variables for the bulk of dimension d. 
We have then argued that this physical interpretation means that the relations between the tails exponent (?],??') and 
the droplet exponent lo should actually remain true on arbitrary lattices when expressed in terms of the dimensions 
(ds , d) , namely 

ds 



ds - 



UJ 



^ (92) 
ds — UJ 

The directed polymer corresponds to the case of a linear object ds = 1 embedded in a space of total dimension d, 
whereas the ferromagnetic random Potts model corresponds to an interface of dimension ds = d — 1 in a space of 
dimension d. These two cases merge for the special case of a linear object ds — 1 embedded in a space of total 
dimension d = 2, which is not surprising since the directed polymer model was precisely invented to model a one- 
dimensional interface in two-dimensional ferromagnetic spin models at low temperature [i^l. This special case also 
explains why it is the distribution of the ground state energy of the directed polymer model (Eq. [T]) which is in direct 
correspondence with the distribution of excitations in ferromagnetic spin models (Eq. [7]). 
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The case of spin-glasses is different for at least two reasons. First of all, only the right tail exponent Vj exists, because 
the distribution of the energy of excitations extends down to i? = as a consequence of the symmetry J — > — J (see 
the discussion around Eqs [75l and [76]) . Secondly, in real space, the dimension dg is expected to be different from the 
value (d — 1) and to reflect the fractal nature of the droplet boundary [1]. 

As a final remark, we should stress that the tails exponents discussed here concern the universal scaling distributions 
of the rescaled variables. But of course, as in the Central Limit theorem, non-universal tails could also be present in 
random systems with particular initial disorder distributions. 



APPENDIX A: TAIL EXPONENTS OF THE GROUND STATE ENERGY DISTRIBUTION FOR THE 

DIRECTED POLYMER ON THE DIAMOND LATTICE 

In this Appendix, we derive the relations between these tail exponents 77(6), 77' (6) defined in EqslMland the fluctuation 
exponent We start from the recursion Eqs fT9l and [20l in the scaling regime. The convolution relation of Eq. [21] 

is simple in Fourier 

GM^\h{q)\ (Al) 

with 

r+°° /■+°° dn 



00 
+ 00 



27r 



r°° da 
Fb{q) = / dxe'9-n,(x) n^ix) = / ^e''"^^ Fhiq) (A2) 

J-oc J-00 

but the recursion relation Eq. [201 is non-local in Fourier. This is why it is difficult to obtain an explicit solution for 
the probability distribution II),. In the following, we show that the problem simplifies for the tails of the distribution 

1. study of the left-tail form Ilb{x — > —00) 

We write the left-tail of 11;, (x) as 

Ub{x) ~ g-c|xr+$(|.|) (^3) 

X — > — 00 

where the function $(|a;|) is subleading with respect to the first term of order This left-tail will determine the 
asymptotic of the Fourier transform Fb{q) of Eq. IA2I for q = is with real s — > -l-oo 

ilb{q = is)= / dxe-'^'Wbix) ~ / da;e"^-'=^''+*(^) (A4) 



s — ^+00 



Since $(a;) is subleading, we perform a saddle point calculation with the two first terms, yielding the saddle value 

x,{s) = {—] ' (A5) 



crj 



One obtains 



y C77(r/ - l)xl (s) 

a. Use of the convolution equation 
The convolution equation is simple in Fourier ( Eq. lAip 

Gb{q^^s)^{Mq = ^s) f =— ^ . 2c(,-i).;(.)+2^(., (^)) (A7) 

- l)xi (s) 
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b. Use of the recursion equation 



We now consider the recursion of Eq. [20] in the tail u — ^ — oo. Using the normahzation 



/+00 r+QO 
dxlli,{x) = 1 = / dxGb{x) 
-oo J —oo 



we may rewrite it as 



1- / dxUbix) 



1 - / dxGbix) 



At leading order in the tail, one has 



The identification 



becomes in Fourier 



dxllb{x) — hi dxGb{x) + .. 



\Tlb C^) ~ bGb{i 



tlb{Q — ^As) bGb{q — is) 



(A8) 



(A9) 



(AlO) 



(All) 



(A12) 



c. Identification 



Rewriting Eq. IA12I using Eq. IA6I and Eq. IA7I yields the following constraints. The identification of the leading 
term in the exponential yields 



and with the notation A = 2" (Eq. (TS] ) this gives 



A~ = 2 



1 -uj 



The identification of subleading terms is compatible with the power-law form 
with the parameters 



1 / crj{r] — 1) 



^ 6V 2^ 
In Fourier, this corresponds to pure exponential forms 



s— ♦+00 o 

s—>--\-oo 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 
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2. Study of the right-tail form Ilb{x +00) 
We write the right-tail of Ilb{x) as 

X — ^+00 

where the function ^{x) is subleading with respect to the first term of order . This right tail will dominate the 
Fourier transform for q ~ —is with real s +00 

/+00 p , 

dxe'^'Whix) ~ / dxe"^-'*^'' +*(^) (A19) 
-00 s^+ooj 

Since ^'(x) is subleading, we perform a saddle point calculation in the two first terms, yielding the saddle value 

"+(^)=('i^V'" (A20) 



One obtains 



flbiq = -t.s) ~ =J (A21) 

dif{if — l)x\ (s) 



s — >+oo 



a. Use of the convolution equation 
The convolution equation is simple in Fourier ( Eq. lAip 

c;j^'(^' _ i)x\ (s) 

The asymptotic behavior of Gb for a; — > 00 is then 

Gb{x) ~ e-D^'''+p(^) (A23) 

with the correspondence 



D =d 21-"' 

p{2x)~2-i'{x) _ 



(A24) 



dr)'{r)' - l)a;'''-2 
This suggests the power-law forms 

e*(") ~ Bx-^ 

^P(x) ^ Qr^'-y (A25) 
with the following relations between exponents and amplitudes 

(T = 2k + — 

6. Use of the recursion equation 

We now consider the recursion of Eq. [20] in the tail u +00. We may use Eq. IA18I and use a saddle-point 
calculation at the left boundary xieft — u/X yielding 

dxllbix) ^ — ^ (A27) 
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Similarly from Eq. IA23I 



+00 



dxGb(x) 



-Du^ +p{u) 



So Eq. [20I becomes 



^-bDu'^ +bp{u) 



(A28) 



(A29) 



c. Identification 



The identification of the leading term gives using Eq. IA24I 

. ^ ^ ^ f 

hD b 



Comparison with Eq. I A13I yields that the two exponents 77 and 77' are related via 

2"^ = (26) 



The subleading terms yield 



Using the power-law forms of Eqs IA25[ one obtains via identification 

K~ba =(fe-l)(77'-l) 



Consistency with Eq. IA26I vields 



1 



2b- 1 
1 

2b- 1 



-l{v'-2)-{b-l){v'-l) 



1 



(r7'-2)-2(6-l)(r;'-l) 



(A30) 



(A31) 



(A32) 



(A33) 



(A34) 



APPENDIX B: REMINDER ON ZHANG ARGUMENT FOR THE LEFT TAIL OF THE DIRECTED 

POLYMER 

Let us now recall Zhang's argument [6j that allows to determine the exponent ?/ of the left tail of the free energy 
distribution of the directed polymer 



(Bl) 



Positive moments of the partition function can be evaluated by the saddle-point method, with a saddle value F* lying 
in the negative tail (jBip 



= y dFPL{F)e-''^^- ~ J dFe-i^y e-'l"^- ~ e<')L^ 



(B2) 



Since these moments of the partition function have to diverge exponentially in L, the exponent 77 of the tail (|Bip 
reads in terms of the droplet exponent 
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This argument can be extended to the free-energy distribution of other random systems. However, as recalled in 
the introduction, the distribution of the free-energy over the samples is simply Gaussian with rj = 2 and iOf = d/2 for 
spin models in any finite dimension d. The distribution of the ground state energy is non Gaussian for the mean-field 
Sherrington-Kirkpatrick model of spin-glasses and we refer to [43| for a discussion of the corresponding tail exponent 
V- 
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